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Abstract: The authors in a previous paper proved the hydrodynamic incompressible limit in d > 3 
for a thermal lattice gas, namely a law of large numbers for the density, velocity field and energy. In this 
paper the equilibrium fluctuations for this model are studied and a central limit theorem is proved for a 
suitable modification of the vector fluctuation field C(0; whose components are the density, velocity and 
energy fluctuations fields. We consider a modified fluctuation field £ e (i) = exp{— e~ 1 tE}( : e , where E is 
the linearized Euler operator around the equilibrium and prove that £ £ (t) converges to a vector generalized 
Ornstein-Uhlenbeck process which is formally solution of the stochastic differential equation d^(t) = 
N^(t)dt + BdWt, with BB* = —2NC, where C is the compressibility matrix, iV is a matrix whose entries 
are second order differential operators and B is a mean zero Gaussian field. The relation —2NC = BB* 
is the fluctuation-dissipation relation. 
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1. Introduction. 



In this paper we study the equilibrium fluctuations for the stochastic lattice gas intro- 
duced in [BEM] . It is a model of particles with discrete velocities jumping on the lattice: a 
particle with a given velocity moves on the 3-d lattice as the asymmetric simple exclusion 
process with the jump intensity chosen so to have a drift equal to its velocity. In each site 
particles collide exchanging velocities in such a way to conserve the number of particles, 
the momentum in each direction and the energy. This model generalizes the one in [EMY2] 
to include the case of particles with different kinetic energy. In [BEM] it has been proved 
the law of large number for this model in the following form. We choose as initial state a 
Bernoulli measure with density, momentum and energy small perturbation (of order e) of 
constant profiles. Then the empirical fields v%(x, t) = isp(e~ 1 x, e~ 2 t), (3 = 0, . . . , 4 of den- 
sity, momentum and energy converge weakly in probability as e goes to to the solution 
of the hydrodynamic equations for this model, which are the incompressible Navier-Stokes 
equation for the velocity field and a diffusive equation (including the transport along the 
velocity field) for the energy. The dissipative terms in these equations are given in terms 
of a diffusion tensor (3, v = 0, . . . , 4, a, 7 = 1, . . . , 3, which is expressed by the Green- 
Kubo formulas. The next natural step is to prove the space-time central limit theorem, 
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namely that the fluctuations fields starting from the equilibrium state converge to a sta- 
tionary multi-dimensional Gaussian process with a given space-time covariance. Since the 
macroscopic behavior of this model is very close to the real hydrodynamics we face the 
main difficulty of the hydrodynamic fluctuations: the Euler terms and the Navier-Stokes 
terms live on different time scales. The same feature is responsible for the impossibility of 
obtaining the compressible dissipative hydrodynamic equations as scaling limit. In fact, 
the previous result on the law of large number is true for an initial condition which is 
a small perturbation of the global equilibrium. This perturbation remains small at later 
times of order e~ 2 and evolves macroscopically according to the incompressible hydro- 
dynamics. The case of the fluctuations is different because a small perturbation of the 
equilibrium may become very large and be of order e _1 on times of order e~ 2 . We go now 
in some details to explain better this point. The fluctuation fields under diffusive scaling 
are defined by 

x 

j3 = 0, • • • , 4, where Gp are suitable test functions, r} t {x) is the configuration in x at time 
t and E is the equilibrium expectation. Ip are the quantities conserved by the dynamics, 
total number of particles, total momentum and total energy in x. 

At time zero the limiting fluctuation fields 

lim$(0,G) = C/j(0,G) 

are jointly Gaussian with covariance 

E[C^(0, G)C„(0, H)] = Cp, v J (fx G{x)H{x). 

The matrix C = (Cp jU ) is called the compressibility matrix. The limit is in the in sense of 
weak convergence of path measures. 

It is not hard to show (it is indeed a by-product of the results and estimates in this 
paper) that the equilibrium fluctuations under Euler time scale are trivial in the sense that 
they satisfy in the limit a deterministic equation. This is a general feature first showed in 
[GP], [FF]. More precisely, the limiting field C E = {(p),P = 0, • • • , 4 

( e (t) = lim Cf (r, G) = lim C e (er, G) 

is solution of the deterministic equation 

d( E (r)=EC E (r)dr, (1.1) 

where E is the linearized Euler operator around the global equilibrium. Equations (1.1) are 
a system of linear hyperbolic equations. The stochastic noise should appear as a correction 
of order e as 

dCf(r) = (A + eV)( E {r)dT + ^BdW T + 0(e 2 ), 
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where V is the linearized Navier-Stokes operator around the global equilibrium and 

BB* = -2VC 

is the fluctuation-dissipation relation. Hence, to see a finite noise one has to look at longer 
times r = e~ 1 t. Formally, since C s (t) = Cf we get 

dC(t) = (e^E + V)C(t)dt + BdW t + 0(e). 

Then the limit lim £ ^oC e (^) does not exist because the Euler modes are too big on this 
time scale. A similar difficulty is present also in the case of ASEP but the analogous of 
E is simply an operator of the form v ■ V x with Vi = (pi — qi)(l — 2a), pi, qi the rates 
of jumping to the left and right respectively and a = E[?y]. Therefore, a Galilean shift is 
sufficient to remove the divergence and in fact in [CLO] the central limit theorem is proved 
for a fluctuation field of the form 

Y e (t, G) = eiJ2 °( £X ~ £ ~ ±vt ) [(Ve-»t(x) - E[rj}] . 

X 

In our case a possible way to subtract the Euler modes is to consider a modified fluctuation 
field which moves together with the waves solutions of (1.1), traveling with velocity of order 
e -1 . Denoting by E* the adjoint operator of E, we define the fluctuation field as 

We prove that the limit e — > exists and satisfies a suitable stochastic differential equation. 
Before writing the equation, we consider the same problem in a very simple case: let A 
and M be K x K matrix with complex entries such that A = —A* where the adjoint is 
relative to the scalar product in M. K . Consider the linear ODE system 

x e = (e~ l A + M)x £ , x s (0) = x. 

Then, y e = e~^ A x e is solution of 

y e = e-i A Mei A y £ , y s (0) = x. 

Consider the limit 

l. 

U := lim I dse~i A Mei A = lime [' ' dse~ sA Me sA 
^°Jo Jo 

An asymptotic average theorem [EP] states that for any 5 > and T > there exists 
£o > such that the solution z of 

z = Uz, z(0) = x 

satisfy 

sup \y £ — z\ < 5, < e < Eq 

0<t<T 
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Therefore, y = lim e ^ Ve is solution of 

V = Uy, y(0) = x. 

The limit U can be characterized in the following way: Let Af be the space of the K x K 
matrices with complex entries. Af is a Hilbert space under the inner product 

(A B) = J2 A li"o- 

l<i,j<K 

For A e Af define Ha as the orthogonal projection onto the subspace of the matrices which 
commute with A 

{BeAf: [B,A] = 0}. 
Since the spectrum of A is imaginary it is easy to see that 

U = lim - f l ' dse~ sA Me sA 
e Jo 

is the projector HaM. 

Applying this kind of considerations to our problem, by Fourier analysis, we prove that 
£(£, G) = lim e ^o £ e (£, G) exists and is a stationary generalized Ornstein-Uhlenbeck process 
characterized formally by the stochastic differential equation 

d£(t) = N£(t)dt + BdW t , 

with BB* = —2NC, where C is the compressibility matrix, N is a second order dif- 
ferential operator and BWt is a mean zero Gaussian field. In particular, this proves the 
fluctuation-dissipation relation —2NC = BB* for this model. Denoting by M the Fourier 
transform of a 5 x 5 matrix whose entries are differential operators, we can characterize N 
as follows: 

N = U & T>, 

Ha the projection on the space of the operators commuting with A. To conclude, we want 
to stress that this procedure of subtracting the Euler modes works in this case because the 
equations for the equilibrium fluctuations are linear. 

The central limit theorem for equilibrium fluctuations is a well investigated topics [S], 
[KL] . A standard procedure is to establish first the tightness of the sequence of fluctuation 
field. Then, the study of the martingale problem allows to identify the unique weak limit as 
a generalized Ornstein-Uhlenbeck process by the use of Holley-Stroock theory. It is crucial 
to evaluate some expression in the martingale problem in terms of the fluctuations field. 
This step, called Boltzmann-Gibbs principle, was first achieved in [BR] for symmetric zero 
range process. The alternative method by [CY] and [CI] was extended to non gradient 
systems by [L] and [C2] . The extension is based on a suitable modification of the fluctuation 
field by adding lower order terms, determined by identifying the diffusion coefficient in the 
hydrodynamic equations. In [CLO] this approach has been extended to a non-symmetric 
case by proving a stronger tightness result and as consequence a stronger Boltzmann-Gibbs 
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theorem. We extend the results on tightness and Boltzmann-Gibbs theorem in [CLO] to 
the present model. Moreover, we prove the convergence of the time averages of the form 
appearing in the martingale problem, by using and adapting some results in [EP] which 
studied the convergence of solutions to the linearized Navier-Stokes equations of solutions 
to the linearized Boltzmann equation. The paper is organized as follows. In Section 2 we 
define the model and recall the previous results on the hydrodynamic limit that we will need 
in the sequel. In Section 3 we define the fluctuation field and state the results. In Section 
4 we identify the limiting distribution of Q e by using Holley-Stroock characterization of 
Ornstein-Uhlenbeck processes with martingales. The Boltzmann-Gibbs principle is proved 
in Section 5 together with the tightness of the process. The theorems stating the existence 
of the time averages are in the Appendix. 

2. Model and hydrodynamic limit. 

We consider the following model introduced in [BEM], which is a generalization of the 
model in [EMY2]: given a finite set of velocities V C M 3 , particles with velocity v G V 
evolve on the sub-lattice Ax, = {— L, . . . , L} 3 , with periodic boundary conditions, according 
to an exclusion process. Collisions between two particles can also occur provided that the 
momentum and the kinetic energy are conserved. The set V is chosen in the following way: 

V = Vi U V 2 , 

where Vi is made of 8 velocities given by 

V = {(±1,±1,±1)} 

and V2 contains 24 velocities, given up to permutation by 

(±w,±l,±l), (2.1) 

where w is some irrational number suitably chosen. 

Formally, if we denote by rj(x,v) G {0, 1} the number of particles on site x G Ax with 
velocity v G V, then the infinitesimal generator of the dynamics on the space Ox = {v = 
(rj(x, v), x G Ax,, v G V)} is defined as 

£ = £ ex + £ c , 

where £ ex is the generator of the nearest neighbor exclusion process with different colors 
(velocities) and C c the generator of the collision process. For a local function / on Ox, 
£ ex is given by 

vev | e |=i xeA L 

where e is a unitary vector of Z 3 (e a , a = 1,2,3, will denote the unitary vectors with 
positive coordinates), r] x ' x+e ' v is the configuration obtained after exchanging the values of 
r](x, v) and rj(x + e, v) and x is some positive constant large enough such that the jump 
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rate is positive. Note that it is chosen so that the drift of the particles with velocity v is 
exactly v. 

The collisions generator C c is given by 

£ c m= E EM^) -/Ml' 

xeA L q eQ 

where Q is the set of admissible collisions, namely the set of velocity quadruples q = 
(v, w, v' , w') G V 4 such that v + w = v' + w' and \v\ 2 + \w\ 2 = \v '\ 2 + \w'\ 2 , and r) x ' q is the 
configuration obtained after the collision on site x between two particles with incoming 
velocities v, w and outgoing velocities v' , w'. Notice that in order to preserve the exclusion 
rule, we take rj x ' q unchanged with respect to r\ if one of the conditions r\[x, v) = 0, n(x, w) = 
0, rj(x, v') = 1 or rj(x, w') = 1 is fulfilled. 

We denote by rj x = {r)(x,v), v G V} the particle configuration in x G Ax,. For a 
configuration rf, the mass, momentum and kinetic energy in site x are 

vev 

I a (Vx) = E( v ' e «)^( x ' v ) ' a = 2 ' 3 ' 
vev 

h(Vx) = E ^\ v \ 2 v(x,v). 
vev 

It is easy to check that the quantities Np(rj) = Yl x Ip(Vx)i P = 0, . . . , 4, are conserved by 
the full dynamics. It is shown in [BEM] that, by choosing suitably the parameter w in 
(2.1), they are the only conserved quantities, in other words this model has the property 
of local ergodicity. 

As a consequence, the grand canonical measures below are invariant for C 

d+l 

VL,n(v) = Z lX II eK p{^2 n (3 I l3(Vx)}, (2-2) 

xeA L 0=o 

where n = (no,..., 714) G M 5 are the chemical potentials and Zl,u is a normalization 
constant. All these product measures are absolutely continuous with respect to the refer- 
ence measure \i obtained by taking n as n := (r, 0,0, 0,0). We set mp = E M [/ / 3(?7o)] for 

(3 = 0, ... ,4 (notice that mp = if j3 = 1, 2, 3) and Ip = Ip — rap. 

In the sequel we call e = L _1 . The law of the process (rj t (x,v)) with generator e~ 2 C 
starting from ji is denoted by and the corresponding expectation by E£. We also call 
fo(v) = E At [r/(a;, v)] the density of particles with velocity v G V with respect to the reference 
measure fx. For any function h on V, we put (/i) = J2 v ev h( v )- 

The currents a of the conserved quantities Ip, (3 = 0, . . . , 4, at site x in direction e Q , 
a = 1, 2, 3, are defined by 

3 

CIp(Vx) = E V a«€, a . 
a=l 
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where, if g is a function on A^, 

Vaj(i) = (V Q £r)(> - e a ) and V a g(x) = g(x + e a ) - g(x). 

Since the local quantities Ip{r} x ) are conserved by the collision generator, we have CIp(j] x ) = 
C ex I/3(r] x ) and the currents can be written as the sum of a symmetric and an antisymmetric 
parts 

< a = xV a / /3 (? ?x )+«;i^ /3 

and 

= (vabxM), # = <^/3^,c»>, 0= 1,2,3 

with 

6x,a(u) = V(x + e a , u)t/(x, u) - + e Q ) + 7/(x, u)). 



Let £ be the space of local functions /ion Al such that 

dW L - n [h] | 



E fi [h]=0 and 



i n=no 



= 0, = 0,. ..,4, 



(2.3) 



dmp(n) 

where mp(n) = W L - n [Ip]. In view of the application of the Boltzmann-Gibbs principle, it 

is important to modify the currents wi^a so that they are in the space Q. It is enough to 
subtract suitable combinations of the conserved quantities and we now get their explicit 
expressions. 

Let n be the chemical potential n = no + 8n = (r + Suq, Sni, Snz, Sn^, 9 + Sn^), then 



1 



l h\)5n a + o(8n), 



W L, n [ w («J.f ] _ C P = 6a/3 [-(\ v \* hl )6n + -<M 4 )Mn J 8n a + o(8n), [3 = 1,2, 3, 



where h = / (1 - / ), ^i = fy)(l - 2/ ) and 

If we denote by Smp = E ML -~ [J/3 (770)] — we get 

<5n = ^((\v\ 4 h )8m - 2(|i>| 2 /io)5m 4 ) , 



(2.4) 



5n 



-8m r 



v\ 2 h ) 

Sn^ =— (2(ho)8rri4 — (\v\ 2 ho)8mo j , 
7 



where 

d> = (M 4 /io)(/i )-(M 2 /io} 2 >o. 



So, defining for ft, v = 0, . . . , 4 and a = 1, 2, 3, 



I n=np ' 



<9m,y(n) 
we obtain 

= M/3,o<W + M/M'W + l{i,2,3}(/3)<W [ao8v,o + 04^,4] , (2.5) 



with 



60 = 3*' 64 = 2 3$' ao= <H^>' a4= <H^>' (2 - 6) 

$ 1 = </i 1 |t;| 4 ></io>-</iik| 2 ></iok| 2 >, 
$2 = </ioM 4 )(/iiM 2 ) - </ii^| 4 )</iok| 2 >. 



Therefore the local function 



= - eg - ^dffK^OjJ + J,^))) (2.7) 



2 ,-0 



belongs to 

Slow-fast modes decomposition of the currents 

We denote by l[ = (io,£, 
over the block of length £: 



We denote by l[ = (1 0,1, ■ ■ ■ , i^) the empirical averages of the conserved quantities 



The measure /i^ m , m 6 I 5 is defined as the canonical Gibbs state of (2£ + l) 3 sites 
with parameters such that l\ = m. It is the uniform probability on the set of 
configurations on the block A^ such that l\ = m. We denote by ag(g) the conditional 
expectation of g given the averages I £ 

a l {g) = W[g\I+]. 

We call C S} £ the symmetric part of the generator £ restricted to the block A^. Since the 
measures /t^ m are the only extremal invariant measures for £ Sj e, we can define £~\g for 
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any function g such that a^{g) = 0. Given any local function g on Q#, the finite volume 
"variance" Vi(g,n) is 

where r. is the translation operator on Q^, r x g(rj) = g(r x r]), t\ = £ — £ 1//9 , £ large enough. 
The "variance" V(G,n) of G is given by 

V(G,n) = limsupV«(G ! ,n). (2.8) 

With an abuse of notation, we denote Vg{G,n) by Vi{G,r,6) when n is the chemical 
potential n = (r, 0, 0, 0, 6>). 

We state here the results in [BEM] 

Theorem 2.1 . There exists a rank 2 tensor D = (-D^) (D^ ,v positive definite matrix) 

and a sequence of local functions = (h^P^ , ol = 1, 2, 3, j3 = 0, . . . , 4) in Q such that, 
setting 

u { ^{ri) = <&) - EE^V^O)) - Z^fa), (2.9) 
where g^(n) is defined in (2.7), it results 

3 4 

limsup^^y^'V,^ = 0. 

1^°° a =l/3=0 



Above Theorem actually holds for any function in Q. 
Lemma 2.2 . The tensor D satisfies 

a-(DC)a= lim E^lTia ■ h^)(-C s )(a ■ h^)]. (2.10) 

q— >oo 

In this formula, a-b = Y^ a =i S/3=o ft a^ ^(^) = Yl x T x9, C-s is the symmetric part of £ in 
L 2 (//), C is the 5x5 compressibility matrix (see (3.3) below for an explicit expression) and 
DC is the tensor {DC)?'" = {D an Cf' v . We define D = D + x l where 1%% = 5 a ,7<W- 



Hydrodynamic limit 

Given functions np(x), (3 = 0, . . . , 4, we consider the Gibbs states with chemical poten- 
tial n(x) = (no(x), . . . , n^x)) 



xeA L (3=0 
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Now, assume that the initial distribution of the particles is ^L,n with n = {rip) the 
slowly varying chemical potentials given by 

n p (x) = \ { ° ] +e\ ( p\ex)+e 2 \f(ex), (2.11) 

where A^ -* = (A^° ) = no and Xp are smooth functions on the 3-d torus T3. We define the 
local equilibrium measure as the Gibbs states HL,n(-,t) with n(-,t) the chemical potential 
given by 

np{x,t) = A^ 0) + sX ( p 1 \ex,t) + e 2 X {2) (ex,t). (2.12) 
Furthermore, we assume 



div A (1) =0, < h x v 2 > + i < h lV 4 > A 



(i) 



0. 



Then in [BEM] (see also [EMY2]) it has been proved that the the law of the process 
at time t > is well approximated by the local equilibrium in the sense that the relative 
entropy per unit volume of the non-equilibrium measure with respect the local equilibrium 
times e~ 2 vanishes in the limit e — > 0. 

We can now state the result proved in [BEM] on the hydrodynamic limit. Let u(z,t), 
z G T3, t G [0, to]) *o > 0, be the classical smooth solutions of the following Navier-Stokes 
equation 

div u = 0, 

A (2-13) 
d t up + dpp + Ku ■ Vup = D a ^d 2 a up, P=l,2, 3, 

a=l 

with initial condition u a (z) = E flL > n <- £ ~ lz '>[I a (rio)] and let £(z,t) be the solution of the 
energy equation 

8 3 

-£ + Hu-V£ = J2 JC «(d 2 a £), (2.14) 

a—l 

with initial condition £{z) = E Mz -'™( £_1 *) [14(7/0)]- The constants appearing in (2.13) and 
(2.14) are given by 



K = 18- 

(h \v\*) 

with h 2 = - 6/o(l - fo)) and 

1 tti - 2C^ 2 
~ (h \v\ 2 ) $2 + ' 

where 



2 ' 



l(hr 


v\ 


4 > 


2 </ii| 


v\ 


2 > 



*i = (h 2 




H 2 >- 


</i 2 


M 4 X*iM 4 > 


#2 = (h 2 


M 4 ><*i 


H 2 >- 


</* 2 


M 2 >(/*iM 4 > 
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Let ¥^ L ' n be the law of the process rjt(x,v) with generator e 2 C starting from the 
measure Pl,u defined in (2.11), with chemical potentials n a (x) of the form (2.12). The 
density (z^j(t, z)), the momentum ((vp(t, z))/3=i,2,3) an d energy (vl(t,z)) empirical fields 
are defined as 

v £ p(z,t)=e 2 6(z-ex)Ip(r] t {x)), 
xeA L 

where Ip(r] x ) = Ip(Vx) ~ mp, mp = E»[Ip(rjo)] and rj t (x) = {n t (x,v), v G V}. 

Theorem 2.3 The density, momentum and energy empirical fields converge, fort < to, 
weakly (in space) in Pg i,n probability, to p(z, t)dz, u(z, t)dz and£(z, t)dz, where ap+bS = c 
for suitable a, 6, c. 

Note that the transport coefficients -D a ,/3 and K, a are suitable combinations of the dif- 
fusion coefficients D@'" in Theorem 2.1. The explicit expressions are given in [BEM], but 
we omit them because they do not play any role in this paper. 



3. Fluctuation field and results. 

In this paper, we are interested in the equilibrium fluctuations of the mass, momentum 
and energy fields. The initial fluctuations, distributed in terms of the measure p, are finite 
but they may become infinite at later very long times because of the effect of waves moving 
with velocity e -1 , which are the solutions of the linearized (around the equilibrium) Euler 
equations (linear hyperbolic equations) for this model. To remove the diverging terms we 
have to modify the usual definition of fluctuation fields not simply by a shift but considering 
fluctuations which move together with the traveling waves 

We denote by Uf the operator exp(— ^E*) where E is the linearized Euler operator, a 
5x5 matrix whose entries are first order differential operators with constant coefficients, 

/ -a Q d \ 
E = -b d -b 4 d 
\ -a 4 a / 

and * is the adjoint with respect the usual scalar product in L 2 (T3,IR 5 ) (the constants a\ 
and bi are defined in (2.6) and d = (<9i, d^-, dz) is the gradient operator). 

For any smooth function G = (G/3)^=o,...,4 : T 3 — > R 5 consider the (scalar) fluctuation 
field £ e on the state space (T 3 )® 5 

T (t, G) = e"/ 2 J2(U!G)p(ex)Ip(r, t (x)). (3.1) 

(3=0 x 

It is equivalent to consider the vector fluctuation field (£|)/3=o,...5 on T 3 whose compo- 
nents (tip) (3=i,... 3 an d ^| are respectively the density, momentum and energy fluctuation 
fields, defined as 

m,<p) = ti e (t,GW), /3 = 0,...,4, 
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where is the vector function with only the j3 component non vanishing and Gp = ip. 

We want to study the evolution of the fluctuation fields in the limit e — » when the 
fields are initially distributed with the equilibrium measure //, given by (2.2). We notice 
that the initial covariance of the limiting fields lim £ ^ £|(0, <p) = £/3(0, (f) is 



^(0,v?)^(0,V)l =C/3,, / (3.2) 



where C is the compressibility matrix (5x5) 



C 



( (ho) (ho^)\ 

%(\v\ 2 h )h 



(3.3) 



with I3 the 3x3 identity matrix, ho defined in the paragraph before (2.4) and ( • ) in the 
paragraph after (2.2). 

Remark that E is not anti-hermitian in L 2 (T 3 , M 5 ), since ao 7^ bo and 7^ 64. However 
a straightforward computation shows that EC satisfies EC + CE* = 0. 

We want to show that the fluctuation field converges to a stationary Gaussian vector 
process with a given covariance. The equal time covariance is exactly (3.2) because of the 
stationarity of the limiting process. 

To state the results we need some extra notation. We introduce the Hilbert spaces 
Tik, fceZ defined by the scalar product 

(G,H) k = (G,L*H) , 

where L = I — A, A the Laplacian operator and (-,-) is the usual inner product of 

L 2 (T 3 ,M 5 ): 

4 r 

(G,H) = J2 dxGp{x)H p {x). (3.4) 



(3=0 



Denote by || • ||fc the norm of TCk and by H-k the dual of TCk with respect to the inner 
product of L 2 (T 3 ,IR 5 ). The fluctuation field (C e (t))t>o is a distribution valued stochastic 
process taking values in the Sobolev space H-k f° r some suitable ko- Its path space is 
D([0, T], Ti-k ) (T > 0), the space of functions with values in TL-k i right continuous with 
left limits, endowed with the uniform (in time) weak (in space) topology. We call Q e 
the law of (£, e (t))t>o when the process is initially distributed according to the equilibrium 
measure fx. It is therefore a probability measure on the space D([0,T],Tl-k o )- 

By analogy with (3.4), we define for local functions g = (go, ... , g^) on Oi and smooth 
functions G = (Go, • • • , G4) on T 3 

4 

(G, g) Q L = e 3 / 2 E Gp{ex)r x g p . (3.5) 

(3=0 xeA L 
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So the fluctuation field (3.1) can be rewritten as 



Z s (t,G) = (U?G,I( Vt (0))) OL . 



We introduce the linearized Navier-Stokes operator V as 



3 



VG = D a „d a d y G. 



(3.6) 



a, 7=1 



Then N = 7Te(T>) is the operator defined as the limit 




T 



dt exp( — E)Vexp(-E)G = tt e {V)G 



and tte(T>)* the adjoint of N with respect to the inner product (•, -) in L 2 (T 3 ,M 5 ). The 
main result of this paper is 

Theorem 3.1 The probability measures (Q e ) converge weakly in D([0, T], 7Y_fc ) to the law 
Q of the stationary generalized Ornstein- Uhlenbeck process £ with mean and covariance 



where (S t )t>o is the semi-group in L 2 (T 3 ,IR 5 ) associated to tv e(T>)* and C the compress- 
ibility matrix. It is formally characterized by the SDE 



One of the main ingredient needed while studying the equilibrium fluctuations is the 
so-called Boltzmann-Gibbs principle which states that the non conserved quantities arising 
in the conservation laws may be replaced by linear combinations of the conserved ones. In 
the context of a non gradient system, the usual statement is not valid and some corrections 
to the fluctuation field have to be introduced (see [C], [Lu]). The situation in the case of 
an asymmetric system is more delicate since the usual Boltzmann-Gibbs estimate is not 
sharp enough and one has to prove a stronger result ([CLO]). We need to generalize such 
a result to the present setup. Indeed we prove the following 




d£(t) = N£(t)dt + BdW t , 



BB* = -2NC. 
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Theorem 3.2 (Boltzmann-Gibbs principle) Assume that h e Q (see (2.3)). Then, 
for any smooth function G : M, + x T 3 1, the following estimate holds 



lim sup E^ - 



<cT\\G\\ 2 V(h;r,e), 



sup £ 3/2 1 / G(s, ex)r x h(rj, s) ds 
)<t<T V Jo V 

where V is the infinite volume variance defined in (2.8). 

4. Limiting distribution of the fluctuation field. 

The theory of Holley-Stroock [HS] characterizes the law Q of the Ornstein-Uhlenbeck 
process £ described in Theorem 3.1 by the following martingale problem: 

M ± (t, G) = £(f , G) - £(0, G) - fds £(s, tt e (V)*G), 

Jo (4.1) 

M 2 (t, G) = (M 1 (t, G)f + 2t (G,Civ E (VyG) 

are martingales under Q. In this section, we will prove that any limit law Q of Q e satisfies 
(4.1). Therefore from the tightness of (Q e ) (see Theorem 5.5 of section 5), it has to 
converge to Q and Theorem 3.1 follows. 

The processes analogous to (4.1) for Q e are 

Mf(f , g) = e (t, g) - e (o, g) - fds e e ( s , ^ E {vyG), 

Jo 

Mf (£, G) = (Mf(t, G)) 2 + 2t (G,Civ E (VyG) 

and we want to show that these processes are martingales up to some error terms which 
vanish as e goes to 0. Given local functions h = (/i a ) a =i,2,3 = (^a)a=i,2,3 ; /3=o,...,4 £ we 
introduce the modified fluctuation field 

3 

C e (t, G, h) = £ e (f, G) - e ^<a Q (t/ t £ G), ^> i , 



a=l 



where (•, -) L was defined in (3.5). Actually we will choose for h the terms of the sequence 

h^ qS> defined in Theorem 2.1, but we will omit the label q for sake of shortness. It is clear 
that the difference between C e (^ G, h) and £ e (t, G) vanishes in L 2 (P^) with e. Moreover, 
it is well known that the following processes are martingales with respect to the usual 
filtration related to the process (rj t (x,v)) 

Mf(t, G, h) = C e (t, G, h) - C e (0, G, h) - f 7i £ (s, G, h) ds, 

t J ° (4.2) 

Mf (t, G, /i) = (Mf(t, G, M) 2 - / 7l(^ G, /i) ds, 

Jo 
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(4.3) 



with 

7l e (t,G,fc) = (a t + £ - 2 £)c e (t,/i,G), 

72 e (t, G, h) = (d t + e- 2 C) (C(t, h, G) 2 ) - 2(C £ (t, h), G) [d t + e" 2 £) (C(t, h), G). 
We first compute the compensator 7f . Let wi a a = (w£"a )/3=o,...,4- Then 

3 

7?(f,G,*) =(d t (U?G)J( Vt (0))) OL -a' 1 J2(da(U!G),jCh a ( Vt )) QL 

a=l 

3 

" £ " 2 E< V -(^ £(? )' XVJMO)) + 4i{rk)) 0tL + Ml G, h), 

a=l 

where, remembering that dtUf = — e~ 1 E*U^ , 

3 

R 5 (t,G, h) = J2(d*(- E *U?)G, h a ( m )) L . 

a=l 

Now, given D a i = (.Df '~)/3,i/=o,...,4) 7 = 1>2, 3, we add and subtract the term 



(D a „d a (UfG),VAvt(0))) , L 

a, 7=1 

in (4.3). Then jf(t, G, h) is equal to 

(dt(UfG),I(i]t(0))) OL + (V*(UfG),I(i] t (0))) OL 

3 3 
a=l 7=1 

+ R 5 (t, G, h) + Riit, G) + R 2 (t, G) + R(t, G), 

with V* the adjoint in L 2 (T3,R 5 ) of the differential operator V in (3.6), c a = (c^) / 3 = o,...,4 
the equilibrium value of the current Wa'* 3 (see (2.4)) and 

3 

R 1 (t, G) = e- 1 £((0 a - £- 1 V Q )(t/ t £ G),xV Q /(r 7t (0))> O I/ , 

a=l 

3 

ifc(f,G) = E {{e- 1 V--d 1 )d a D a „{u;G), l( Vt (0))) OL , 

a, 7=1 
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3 

R(t,G) =e" 1 ^((da-e-^a) (U t <G)Ai-**)o, L - 

a=l 

From the definition of UfG : the first term of the sum (4.4) can be written as 

3 3 
a=l a=l 

where the coefficients of the matrix d a = (g^' i/ )/3,i,=o...,4 were defined in (2.5). Recalling 
the definition of the local functions g a = (^)/3=o...,4 and u a = ("Uq)/3=o...,4 in (2.7), (2.9) 
(we omit the label q), we obtain 

4 6 

7l e (t, G, h) = (V*(UfG), U Vt (0))) Q L + Ri& G)+J2 Ri(t, G, h), 

i=l i=5 

where 

3 1 
R 3 (t,G) = e- 1 ^((da-e-'Va + -V" d a ) (UfG), dJ(r k (0))) Q ^ 
<■■ i 

3 

R 4 (t,G)=e- 1 J2((d« -e" 1 V a )W,ffa(>? ( )) 0iL 

a=l 

and 

3 

Re(t, G, h) = -e- 1 ^(«9 Q (£/ t e G), u a ( Vt )) Q L . 

a=l 

To prove that the compensator j^{(t,G,h)ds is converging, we have to control the re- 
mainder terms. 

The remainders for i = 1, 2, 3 are easily controlled by the following 

Lemma 4.1 Let h be a mean zero local function and G : M + x T 3 — > R a continuous 
function. Then there exist a constant c depending only on h such that, for all t > and 
a// e > 

<cT 2 ||G||L. 



The proof is an easy consequence of the Schwartz inequality and the stationarity of P^. 
We refer to Lemma 4.1 in [CLO] for details. 

By Taylor expanding and using Lemma 4.1 we immediately obtain 

lim E^ 



sup ( I dse 3/2 ^2G(s,ex)T x h(r] s ) 

0<t<T\Jo „ 



ds Ri(s, G) 
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for % = 1, 2, 3. 

The other terms are estimated by using the refined Boltzmann-Gibbs principle (Theorem 
3.2) because the functions g@ and are in Q (h@ E Q by hypothesis). We get 



lim E^ 

e^O e 



t - 2" 

cis Ra(s,G) 







= lim E^ 

e^O e 



t 

(is -Rs(s, G, fo) 



o 



= 



and 



lim sup E^ 



Of 



(is i?6(s, G, /i) 



3 4 



< 



ct max 3 ||«9 Q (tTG) ||g J] E r ' 

a=l/3=0 



From the definition of the semi-group (U£), it is clear that ||<9 Q (£7. e G)||o is uniformly 
bounded in e. Moreover the diffusion coefficients are chosen in such a way that, 

since we take for hP a the terms of the sequence (ha' 13 ) given in Theorem 2.1, we have 

lim V(u^' ;r,9) = 0. 

We have shown so far that there exists a random variable R% which converges to in 
L 2 (P^) as e — > and then q — > oo such that 

Mf(f, G, /i) = r (f, G) - r (0, G) - fds (V*(U s s G)J( Vs (0))) L + R% (4-5) 

Jo 

We would like to have instead of the third term in (4.5) a term of the form 

(U!(HG),I( Vs (0))) o , L 

for some suitable operator H that we could then rewrite as £ e (s, HG), so to identify the 
limiting martingale problem. We proceed in the following way: 

(V(u;g)J(v.(o))) , l = (ukud-^uigJmo)))^. 

By Lemma A. 2 and Lemma A. 3 



lim E^ 



ds 



e (s, (UI^VUZG) - e {s, n E * (V*)G) 



= 0. 



Hence, noticing that 7Te*(2}*) = 7Te(P)*, we have proved that there exists a random 
variable G| which converges to in L 2 (P^) when e — > and then (/ — > oo such that 

M!(t, G, /i) = e(t, G) - £ e (0, G) - /'ds r (s, tt e {V)*G) + C%. 

Jo 
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We now compute the compensator 7§ in (4.2). We first remark that 

-f e 2 (t,G,h) = (e- 2 C)(C(t,G,h) 2 )-2C(t,G,h)(e- 2 C)( £ (t,G,h). 

We introduce the operator = £ ex '( 2 ) + £ c '( 2 ) for local functions / and g as 

C^ 2 \f-g) = C ex (fg) - fC m g - gC ex f, 
C^ 2 \f;g) = C c (fg)-fC c g-gC c f. 

Then we obtain 

72 e (t, G, h) = Yffa G) + Kf (t, G, h) + Yi(t, G, h), 

where 

4 

Y 1 %t,G) = e E(^ G )/3(^)(^ G )^^) £eX ' (2) (^fe);t(%)), 

/3,i/=0 x,y 

3 4 

Yi(t,G,h) = -2e 2 J2 E E(^ G )/3^)^(^ G )^^) £eX ' (2) (^(^);^a)), 

a=l /3,i/=0 x,y 
3 4 

a, 7=1 /3,i/=0 a;,y 

From the explicit formulas (4.7), (4.8) and (4.9) that we will get below for Y{(t,G), 
1, 2, 3 and the use of Lemma 4.1 it is easy to see that 

'da ( 1 s 2 (s,h,G)-E^[ 1 s 2 (s,h,G)]) 

is converging to in L 2 (P^). So, all we need to compute is [7! (t, ft, GQ] . 
Notice that 

£ ex ' (2) (/;#) = 6(x,x + e, f)V :E , x+e ,„/V :E , x+e) „5r, 

x,e,v 



x,q 



with V x , x+e , v f = f{rf - fir)), V x , q f = f{rf«) - f( V ) and 

b(x, y, v)= (x + \v-{y- x))t/(x, u)(1 - u)). 
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So, if we let 4>o(v) = 1, (f>p(v) = vp for (3 = 1,2,3 and 04 (u) = |M 2 , a straightforward 
computation leads to the following 

3 4 



I X 



l?(f,G)=e 3 £ ^ 52(d a {U!G)pd a (U t e G)v){exy. 

a=ip,u=0 X 

E [&(#, x + e Q , u) + 6 (a; + e Q , x, u) <j)p(v)(j) u (v ) + 0(e). 



(4.7) 



Therefore 



3 4 



a=l P,v=0 x 

= -2 X ((U?G),AC(U?G)) + O(e), 

where A is the vectorial Laplacian operator defined as (AG)p = AG p. Observe that 

((UfG),AC(U?G)) Q = (G,e^ E Ce-^AG) = (G, Ce^ E * e~^ E * AG) Q , 

where we have used that EC = —CE*. In conclusion, Yf(t,G) converges in L 2 (P^) to 
-2 X (G,CAG) Q . 

We get in the same way 

3 4 



Yj(t,G,h)=2e* E 52(d a (U;G)MUZG)u)(ex)x 

a, 7=1 p,v=0 x 

^2 [ b ( x i x + e a ,v)-b(x + e a , x, u)] <i>p(v)V x , x + ea , v T{h v a ) + 0(e), 

V 

(4.8) 

where T(/i^) = J2 x T x h^- Since fx is invariant for the jump generator of particles with a 
given velocity, it is easy to check that 



= 0, 



r](x, v)(l - r](x + e a ,v))V x , x+ea , v T(h u a ) 

which implies that the time integral of Y^t, G, h) converges to in L 2 (P^) by Lemma 4.1. 
The last term Y" 3 £ (t, G, h) is given by 



3 4 



£ " E E 52(d a {u;G) s d^u;G)v){ex)x 



a, 7=1 P,v=0 x 



bfa * + e, v)V x , x+e , v T(hi)V x , x+e , v r(h») + £ V x , q T(hi)V x , q T(h») 



v;\e\ = l 



(4.9) 
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By using again Lemma 4.1 it is immediate to show that the time integral of Y" 3 e converges 
in L 2 (P^) to its average that we are going to compute. 

Let C e s x be respectively the symmetric part of £ ex in L 2 (pt). It is easy to check that for 
any local function / and g 



v;\e\ = l 



and 



»/(0,t;)(l-»7(e,t;))Vo >e> t,r(/)Vo > e > t,r(<y) = 2E^ T(f)(-C?)g 



52w[Vo tq T(f)V 0tq T(gj\ = 2E^T(f)(-C c )g . 



Therefore 



E^ [Yi (t, G)} = 2e 3 ^ [ F i d i U t G ) M • h) (-C 3 ) (d(U?G) (ex) ■ hj\ , 



where C s and a ■ b were defined after (2.10). Remember that the functions h = (h^) are 
chosen as the terms of the sequence (ha ^) in Theorem 2.1. Lemma 2.2 asserts that 



lim W 

q^oo 



T(a-h^)(-C s )(a-h^)\ = 2a-(DC)a. 



Hence, 



W[Yi(t,G,h)} = 2e 3 J29(U t e G) ■ (DC)d(U t £ G) + o q (l) 

X 

= -2(U?G, (VC)(U?G)) Q + o,(l) + 0(e), 
where, denoting by D a i the matrix (-D^)^^...^, 



VG = D an d a djG. 

a, 7=1 

With the property PC = CP*, we get 

(U t s G, (VC)(U!G)) Q = (e-^G, (CV*)e^ E * G) Q = <C, Ce^ E * V*e^ E * G) Q 
and by Lemma A. 2 



lim / 



lim / ds exp(-£*)P*exp(-j£*) = t7r E *(P*) =tn E (Vy 



so 



lim 



dsY 3 e (s,G,h) = -2t{G,Cn E (VyG) Q . 
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To summarize, we have proved that there exists a random variable Rj vanishing in L 2 (P^) 
in the limits e — * and then q — > oo such that 

M|(t, G, h) = (Mf(t, G, h)) 2 + 2t(G,Cn E (VyG) + 2t(G,CAG) Q + R q e 
= (M[(t, G, h)) 2 + 2t(G,Cn E (VyG) + flf . 

This completes the proof of Theorem 3.1, once the Boltzmann-Gibbs principle and Lemmas 
A. 2 and A. 3 are proved. 



5. The Boltzmann-Gibbs principle. 



Since we closely follow the strategy proposed in [CLO] to prove Theorem 3.2, we will 
only focus our attention to the points where non trivial changes are necessary. 

One of the ingredients in the proof is the equivalence of ensembles, which is classical for 
Bernoulli product measures but, as far as we know, is not in our case. We state below a 
weaker statement which will suffice to our purpose. 

For a given chemical potential n G R 5 , let M(n) = (Mo(n), . . . , M^(n)) be defined as 
M p (n) = W L < n [Ip(r]Q)}. If we put A = M(R 5 ), it is easy to verify that n i-> M(n) is a C 1 
diffeomorphism from IR 5 onto A, in particular the inverse function M i— > n(M) is continuous 
on A. Given a > 0, we introduce the set A a of M G A such that, |n(M) — uq\ < a, with 
no = (r, 0, 0, 0, 0) the equilibrium chemical potentials. We denote by J1l,m the grand 
canonical measure Hl,u(m) which satisfies therefore E^ L ' M [Ip(r]o)] — Mp for f3 = 0, . . . , 4. 

Recall that I^iv) = (/o ■ ■ ■ > it hi 1 !)) are empirical averages of the conserved 

quantities in A/,. For any particle configuration r\ in O^, we call N^(rj), v G V, the average 
number of particles with velocity v in A^. 

Also recall the definition of <j>p(v) before (4.6). Given k = (k v ) v£ \>, we set Ip(k) = 
^ v </> p (v)k v ™dl+(k) = (l (k),...,h(k)). 

Lemma 5.1 (Equivalence of ensembles) Let h be a local function. Then there exists 
a constant c = c(h, a) such that 



E»[h\l+ = M] -E^ L - M [h] 



< ce- 



uniformly in M G A a . 

Proof. Let t be the number of velocities in V and denote by v a , a = (a„)„ e v 5 the product 
measure on Ql of Bernoulli measures with parameters a = (a v ) v , i.e. E Ua [rj(x, v)] = a v . 
A straightforward extension of the classical strong equivalence of ensembles asserts that 
for any local function g, 



E u «[g\N£ = k v1 v GV] - E Uk [h] < C(h)e 



(5.1) 



uniformly in k = (k v ) ve \; G Bl = {0, L 3 , . . . , 1}^ 
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We first compute the term E M [/i|/^ = M]. Since this expectation does not depend on 

the chemical potential (here no), it is equal to E" 1 / 2 [h\I^ = M] with the obvious abuse 
of notation 1/2 = (1/2, . . . , 1/2). Therefore, from (5.1), 



W[h\l+ = M] 



Z k eB L Jt(k)=M"i/2(N v L = k v ,ve V)E»"[h] 



E 



:w 

keB L ,I+(k)=M 



vi/2(NZ = k Vt veV) 



+ 0(e 3 ). (5.2) 



Since the particles with different velocities are independent 

v 1/2 (Nl = k v ,veV) = H u 1/2 (N v L = k v ) 

vev 

and the asymptotics of a single term in the product above is given by the Stirling formula 



vi/2{N v L = k) 



^/2TTE- 3 k(l - k) 



exp 



-e-3( sW+ l 0g2 )]( 1+0 (^)), (5.3) 



where s(k) = klogk + (1 — k) log(l — k) is the entropy. In particular, if (k v ) v belongs to 
B S L := B L n [5, 1 - S] £ for some small 6 > 0, then 



vi/2(NZ = k v ,veV) = 



1 



y/(27ve-^YYl v k v (l-k v ) 



exp [s- 3 J2(s(k v )+log 2)] (l+0(£ 3 )) . 



The fact that the entropy is convex suggests to use the Laplace method to derive the 
asymptotics of both terms in the ratio (5.2). This is the aim of Lemma 5.2 below which 
is stated in the £ = 1 case without any constraint on k, nevertheless the generalization to 
higher dimension with constrains is easy because, up to a linear change of variables k i— > k 1 ', 
the sums over k in (5.2) with constraints can be written as a sum without constraint over k' 
in a cube of dimension I — 5 (5 is the number of linear conditions I + (k) = M). Therefore, 
we have 



keB 5 I+(k)=M 



E "i/2(NZ = k v ,veV)E'*[h] 

vm^mw^m exp ^ {S{K) + log 2) ] ^ [h] (1 + ° (£3)) ; 



where k* is the 
I+(k) = M, 



of J2 v (s(k v ) + log 2) under the constraints k G [5, 1 — 5] n and 



T P = 



E 



exp 



^ vS " {K) e- 3 (k-k*) 2 



keB s L ,i+(k)=M 

\k-k*\<e' ia 

with < a < 1/2. Notice that this result holds provided that the minimizer k* satisfies 
k* e]5, 1 — S[ l , that will be shown below. As a consequence, the ratio (5.2) is equal to 

E^[h](l + 0(e 3 )), 
22 



provided that the contributions from "bad" configurations are negligible. 

Let k be the minimizer of J2 v (s(k v )+\og2) under the constraints k G [0, 1] , I + (k) = M. 
From Lagrange optimization theorem, k has to minimize the function 

4 

+log2) + J2lpJ2(Mv)k v - M p ), 

v /3=0 v 

where 4>p have been defined in the line before (4.7) and (7^) are Lagrangian multipliers. 
So the minimizer satisfies 

4 

s '( K v) = 7/30/3^), v G V. 
/?=o 

Since the derivative of the entropy s'(a) is equal to the associated chemical potential 
A = log jz^, we have v K = 7 = (70, • • • ,74) but the constraint I + (k) = M implies 

that 7 = n(M) that is to say v K = Jil,m and in particular E UK [h] = E^ L ' M [h]. Moreover, 
if A„ = log 1 is the chemical potential related to k v , then we have 

4 

(3=0 

From the assumption M G A a , the previous equality implies that we can choose 5 > 
small enough such that k e]25, 1 — 25[ e uniformly in M G A a . Such a choice of 5 implies 
that k* = k. 

So the lemma will be proved if we finally show that the contribution of the "bad" k 
(k G B L \ [6, 1 - 8} e ) inside the sums in the numerator and denominator of the ration (5.2) 
is irrelevant with respect to the leading term. From Stirling formula (5.3), there is c > 
such that 

J2 vi/2{NZ = k v ,veV) <exp[- £ - 3 ^( S (/^)+log2)-clog£]. 

k<£B s L ,I+(k)=M v 

From the discussion above, there exists b > such that J2 V s (kv) > J2 V s ( K v) + b, therefore 
"i/2{NZ = k v ,veV) < cexp[- £ - 3 ^(s(^) +log2)]exp[-6 £ - 3 /2]. 

k£B s L ,I+(k)=M v 

□ 

Lemma 5.2 Let ip and cj> be smooth functions on [0, 1], tp concave, non negative. Assume 
that the maximizer 9 of ^ is in ]0, 1[, then 

i=0 
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with 9 the maximizer ofip and 



ex p 

\i — Nd\<N 1 - a 



(i - NO)' 
N 



< a < -. 



Proof. We start by factorizing the leading term exp[Nip(9)] in the sum. For simplicity call 



MO = 0(^)exptfM^)-lK0)]. 



-AT- 



From the assumption on 9, if 5 > there exists a constant c(S) > such that 



V U N (i)<exp[-Nc(S)]. 



i-N8\>NS 



Moreover, choosing 5 small enough ensures that ip(x) — ip(9) < —c(x — 9) 2 provided that 
\x — 9\ < 5, where c > is a constant which will change from line to line. Then, given 

< a < |, 

£ C/ivW^expi-AT 1 - 2 ^)]. (5.4) 

N 1 ~ a <\i-N9\<NS 

So the main contribution is coming from ^| i _ A r 6 »| <A ri-a Un{i). Using Taylor expansion, 
we see that in this range of i's, 



Un(i) = exp 



«o + aA^) + a obl N^)\o((^ 2 ) + o(N(^ 4 



where j = ± -6, a = (f)(9), a± = (f>'(9), b - 2 
oddness" of j, 



< and b 1 = ^H. 



By the "almost 



ex P 

-N6\<N 1 - a 



bo 



N 



ai (i) +aoblN (j^y) =O{^)S N (a,b ). (5.5) 



We also remark 



\i-Nd\<N 1 - a 



; \ 2 



exp 



6o : 



\ J ~ AT 



ex P 



&0 



2N 



so that 



£ (i) ; 



exp 



AT 



= 0(^)Ma,£), 



(5.6) 
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We get in the same way 



\i-N6\<N 1 - a 



exp 



bo 



0{±)S N (a,±). 



(5.7) 



Finally, comparing N 1 / 2 Sn{ol,o) (a < 0) with the integral of a Gaussian, it is easy to 
check that Sn{oc, a) = c(a)VN + 0(1), therefore 



£W (a, y) = 0(l)SW(a!, 6 ). 



(5.8) 



Putting together formulas (5.4) to (5.8), the Lemma is proved. 
□ 

Even if the equivalence of ensembles that we stated in Lemma 5.1 is weaker than the 
classical one, it is enough to prove the following result which is actually the only estimate 
needed in the proof of Boltzmann-Gibbs principle. 

Corollary 5.3 If h e Q is a local function then 



W 



{w[h\it = M]y 



Proof. Let h = E^[h\l^], h(M) = E^ L > M [h] and consider the decomposition 



W[h 2 ] < 2E^(h-~h(I+)) 2 j +2E^[(/i(I+)) 2 ]. 



Since h is in Q, we have 



h(m) = and 



dh 



M=m 



with mp the equilibrium values of Ip{r]o). Therefore 



\~h(lt)\<c \{hAv)-rnp)(h, L {rj)-m„)\ 

0,v=O 

Hence E»[(~h(I+)) 2 ] < ce 6 . 

On the other hand, it results from Lemma 5.1 that for any a > 



(h-~h(I+)) 2 <c(a)e 6 + cP»[l+(£A a ]. 
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From the continuity of the function M i— > n(M), there exists b > such that 



Finally, since l£ = J2 X I + (i] x ) with I + (n x ) i.i.d. random vectors with finite exponen- 
tial moments and expectation m under fj,, a large deviation estimate provides 

P^[\l£ - m\ > b] < exp(-ce 3 ). 

□ 

The first result used in [CLO] (Lemma 4.3) before establishing Boltzmann-Gibbs is a 
general estimate bounding the equilibrium expectation of the squared time integral of zero 
mean functions of Markov processes by their || • norm. More precisely, if X is a Markov 
process on the finite state space S with generator L and ergodic invariant measure tt, then 
there exists a (universal) constant c > such that for any function / : [0, T] x £ — > K. 
satisfying W[f(t, X t )] = for any t E [0, T], we have 



sup 

0<t<T \Jo 



t x 2- 

dsf(s,X s ) 



< c 



ds\\f(s, 



(5.9) 



where 



. 1 =sup{(f,g) + (f,L a g)} 



and L s is the symmetric part of L in L 2 (n). 

The next lemma (Lemma 4.4 in [CLO]) is needed to control remainder terms in the 
proof of the Boltzmann-Gibbs principle. 

Lemma 5.4 For any local function h e Q, there exists a constant c(h) > such that for 
any subset B of Al, any smooth function G : [0, T] x T 3 — > R and e small enough, 



<c(h)(l + T) / dse 3 ^G 2 (s,£x). 

(5.10) 

Proof. Following [EMY2] (section 4), we introduce an alternative representation for the 
particle configuration rj x = (rf(x,v)') v & t site x: one can find n(n — 5) numbers cp(v), 
(3 = — n + 5, . . . , — 1 and v G V, such that if we put 

vev 

then the map r\ x \— > (i/3(f7ai))/?=-n-i-5,...,4 is one to one. Moreover the coefficients cp(v) can 
be chosen in such a way that the covariances K^[Ip(rj x ); I v (r) x )\, (3 ^ v, vanish (except 

[3, v e {0,4}). We also introduce the variables Ip(r] x ) = Ip(r] x ) — E^/g^)]. 
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E^ sup ( [ dse 3 / 2 ' 1 ^ G(s,ex))T x h(r] s )ds) 
o<t<r\Jo ^ B J 



Let Q ex be the space of functions h such that E^[h] = and ^ x E M [/i; Ip(r) x )] = for 
any (3 = —n + 5, ... ,4. The integration by parts lemma valid for ASEP (Lemma 6.1 in 
[EMY1]) easily generalizes to a superposition of ASEP. 

We now turn to the proof. Fix h G Q, we can find coefficients (a ( a) ( a<o such that 
h - Y.p<o a ph,L is in G ex , where Ip^iv) = |Al| _1 Ex Ipiflx)- Therefore, it is enough 

to prove the lemma in the case where h £ Q ex and in the case where h = Ip^iv)- The 
first case is a straightforward generalization of Lemma 4.4 in [CLO] since the integration 

by parts formula is valid in Q ex . In the second case, denote by Ip(i] x ) (/? < fixed) the 
conditional expectation of Ip(rj x ) with respect to the empirical averages of the conserved 
quantities 1^. Then the left hand side of (5.10) is bounded above (up to a factor 2) by 
the sum of the two terms 



sup 

0<t<T 



(5.11) 



where Gf := £ 3 Exes ^( s ' £x )i an< ^ 



sup 

0<t<T 



(5.12) 



From the inequality (5.9), (5.11) is less or equal to 



cV L fy Vo ),r,6) [ ds(Gf) 2 
Jo 



and by corollary 4.6 of [EMY2], V(lp(r]o), r, 9) = limsup L V L (J^(r/ ), r, 0) < +oo. So 
(5.11) is bounded above by 

c dse s y2 G 2 (s,ex) 



(c a positive constant). Finally, by stationarity of \i 1 the term (5.12) is less than 



e-*W[{ip{ m )f]T [ ds{Gf) 2 . 
Jo 



From Corollary 5.3, E^ [(/^(ryo)) 2 ] < ce 6 and (5.12) is going to zero as e — > 0. 
□ 



Finally, Corollary 5.3 and Lemma 5.4 allows to extend straightforwardly the proof of 
the Boltzmann-Gibbs principle given in section 4 of [CLO] and then to obtain Theorem 
3.2. 
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We conclude this section by pointing out that the arguments for the proof of tightness 
(section 5 of [CLO]) can be easily adapted to our case. Notice that, up to now, we did 
not need to have the supremum over time inside the expectation in the Boltzmann-Gibbs 
statement, however it is used in this part to control some terms arising in martingale 
compensators. So, we can state 

Theorem 5.5 The family of probability (Q e ) £ >o on D([0, T], Ti-k ) is tight since 

lim limW sup Ul |U > M) = (5.13) 

M^oo e^O Vn</<T ' 



and for any a > 



limhmW sup ||£ t £ - £f||_ fco > a) = 0. (5.14) 



S^ue^u ^\ s -t\<6 

0<s,t<T 



Appendix 

Let A be the space ofnxn matrices with complex entries. A is a Hilbert space under 
the scalar product 

(X,Y)= Xk,tY k ,i. 

l<k,£<n 

Given a matrix A in A the projector Ha is defined as the orthogonal projection onto C(A), 
the commutator space of A 

C(A) = {M e A: [M,A] = 0}, [M, A] := MA - AM. 

Lemma A.l Let A be a diagonalizable matrix, Sp(A) G iR. Then, for any matrix M and 
t > 

I ft s s 

lim - ds exp (-A)Mexp( — A) = U A (M). 
e^o t Jq e e 

Proof. We follow the proof in [EP] . Let P be a non-singular matrix and R a real diagonal 
matrix such that A = P~ 1 iRP. Let {Sj,j = l,...,m} be a partition of the integers 
{j = 1, . . . , n} such that 

Rk = Ri if k, £ e Sj for some j, 
Rk 7^ Ri otherwise, 

where Rj,j = l,---,n are the eigenvalues of R. We define the bar operation in the 
following way: let K = (Kk,e) € A be 

_ f 1 if k,£ E Sj for some j, 
1 otherwise. 
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Then M , M e A, is denned as 

m m = jr M M M . 

Observe that M is the diagonal part of M in the simple case of R with distinct eigenvalues. 
We have that 

exp(^A)Mexp(-^A) = P' 1 exp( jR)PMP~ 1 exp(- y-R)-P- 
It is proved in [EP] that 

P~ X PMP- X P 

is a projection onto C(A). Hence M = IIr(M) because R is diagonal. Moreover, IIr(M) = 
Ha(M) because R is diagonal. So it is enough to prove that for any M 

rt 



1 f l i s i s 

lim - ds exp(— R)Mexp( R) = M. 

s^o t Jq e e 



In [EP] it is also shown that for any matrix M there exists a matrix S such that M can 
be decomposed as 

M = M+ [S,R]. 

Since M commutes with R 

exp(jR)Mexp(-jR) =M + exp(jR)[S,R]exp(-jR). 
The second term on the r.h.s gives 

(exp(-i?) [S, R] exp(--R)) = S k , t {R t - R k ) exp(-(R k - R e )), 
where S = (Sk,e) and R = (Rk,e) = (Rk$k,e)- As a consequence, 

1 1 S 7 S 

lim - / dsexp(-R) [S, R] exp( R) = 0. 

e^o t Jq e e 



Lemma A. 2 Let E 6e a first order differential operator such that its Fourier transform 
E(fc) satisfies Sp(E(fc)) G zlR for any k and let D = : D cen d ce d 1 be a second order 

differential operator, where D = (D an ) = (D^) is a definite positive rank 2 tensor. Then 
there exists a definite positive second order differential operator 7Te(D) such that for any G 
smooth 



lim 



ds 



exp (-E)Dexp( — E) — 7Te( O) 



= 0. 
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Proof. Let D(/c) be the Fourier transform of D 



a, 7=1 



It is enough to prove that for any t > and for any G smooth there exist a matrix 7Te(B) 
such that 



lim 



(is 



exp (-E)Dexp( — E) -tte(D) 



G 



= 



where || • ||o is the usual norm in L (T 3 ,IR 5 ). Choosing 7Te(D) = 7r t(®)' that * s an eas y 

consequence of Lemma A.l via dominated convergence theorem since, by assumption, E is 
diagonalizable with pure complex eigenvalues which implies ||exp(|E) || < const. Finally, 

since D is definite positive, the same is true for ^(D). 
□ 

Notice that Lemma A. 2 implies that for any < s < t, 



lim 



du 



11 u 

exp (-E)Dexp( — E) - tt e (D) 



G 



= 



Lemma A. 3 Let A s (s), A be linear operators from Hk +2 to Tik such that 

< oo 



sup ||A e (s) 

e,0<s<t 



I fco+2— >fc 



and for any G e Ti-k+2 and < s < t 



lim 



du [A £ (u) - A]G 



= 0. 



fc 



Then, for any G G 7~Ck +2 



dse(s,[A £ (s)-A]G) 



0, 



where £J is the fluctuation field. 

Proof. We set (£ t e , G) = £ e (t, G). Let = t < t x < . . . < t e = t be a subdivision of the 
interval [0, t] of size 5 > 0. Then 



£-1 



ti+1 



[ dse(s,A £ (s)G) = J2(^ I dsA e (s)G)+Rl 
Jo t=0 x J ti 
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with 



Since 



i=0 



\R\\ < t sup ||^ - CJUo sup \\A s (s)\\ ko+2 ^ k JG\\ ko+2 , 

\s 1 -s 2 \<& ' 0<s<t 

0<si ,S2 <t 

it results from tightness (5.14) that for any 5 > 0, 

lim lim P^O^fl > 5) = 0. 

(5^0 e^O ^ v 1 y 



Moreover 



with 



/ <h (C, A s (s)G) = J2(^ A G)(ti+i ~ U) +Ri+ R £ 2 , 



1-^2 I — 



£ f tl+1 ds(e u (A%s)-A)G) 

1+1 <is[A e (s) - A]G 



< t SUp IKfll-fco max 

0<s<t * 



fc 



From assumption 



lim max 

e^O i 



1+1 ds [A £ (s)-A]G 



= 0. 



So, using tightness (5.13), we get for M > and £ small enough 

P £ {\W 2 \ > 5) < P £ ( sup O_ fc0 > A) 

which vanishes in the limit M — > after e — > 0. With the same kind of arguments (using 
tightness again), we get 



where 



J2(^AG)(t l+1 -U)= ds (£, AG) + iC 
lim limPf (imi > 5) = 0. 



We have proved so far that J* * <is(^ e (s, A e (s)G) — £ e (s, AG)) converges to in proba- 
bility. To assert that the convergence occurs in L 2 (P M ) it suffices e.g. to check that 



supE^ 



t ds(^ s ,(A%s)-A)G)\ 



< oo, 



which is clear from the assumptions on the operator A s (s) and 
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